Search for: All records

Abstract When numerically solving Einstein’s equations for the evolution of binary black holes, physical imperfections in the initial data manifest as a transient, high-frequency pulse of ‘junk radiation.’ This unphysical signal must be removed before the waveform can be used. Improvements in the efficiency of numerical simulations now allow waveform catalogs containing thousands of waveforms to be produced. Thus, an automated procedure for identifying junk radiation is required. To this end, we present a new algorithm based on the empirical mode decomposition (EMD) from the Hilbert–Huang transform. This approach allows us to isolate and measure the high-frequency oscillations present in the measured irreducible masses of the black holes. The decay of these oscillations allows us to estimate the time from which the junk radiation can be ignored. To make this procedure more precise, we propose three distinct threshold criteria that specify how small the contribution of junk radiation has to be before it can be considered negligible. We apply this algorithm to 3403 BBH simulations from the Simulating eXtreme Spacetime catalog to find appropriate values for the thresholds in the three criteria. We find that this approach yields reliable decay time estimates, i.e. when to consider the simulation physical, for $>$98.5% of the simulations studied. This demonstrates the efficacy of the EMD as a suitable tool to automatically isolate and characterize junk radiation in the simulation of binary black hole systems. 

ABSTRACT Biomimetic designs are inspired by the complex and unique behavior of naturally occurring materials, and can be applied to many systems, including polymers. ZIPer polymers (Zwitter arene‐ion like polymer) are inspired by byssal threads found on mussels, and their physical state is highly sensitive to various environmental conditions. Specifically, the ZIPer polymer undergoes chemospecific phase transitions, exhibiting potential for its use as an ionic responsive technology. Though this phenomenon has been observed with Raman spectroscopy, little is known about how salt identity or concentration affect polymer inter‐ and intra‐chain interactions. Previous studies have used Raman spectroscopy to analyze ZIPer polymer behavior in the presence of salt; however, the effect is typically only observed with sodium chloride and often only compares spectra at two concentrations. Additionally, studies have mainly focused on the spectral evidence of cation–π interactions, significantly narrowing their spectral range. In order to develop a more predictive framework for ZIPer polymer behavior, a range of salt identities and concentrations need to be tested. This study uses Raman spectroscopy to investigate ZIPer polymer behavior in the presence of a series of salts, namely NaCl, NaOTFA, NaBr, NaBF₄, and NaPF₆, each at 0.1 M, 0.5 M, 1.0 M, and 1.5 M concentrations. Moreover, we observe spectral changes in a range from 550 to 2000 cm⁻¹. Spectral evidence suggests that the cation–π interactions previously hypothesized to be the driver of ZIPer polymer behavior are not the only mechanism determining the chemoresponsive phase transitions. We hypothesize that cation–π interactions and dispersion forces are competing mechanisms controlling ZIPer polymer behavior. Furthermore, we suggest that at certain concentrations the dominating mechanism transitions, and this inflection point is salt identity dependent. 

Abstract Gravitational memory effects and the BMS freedoms exhibited at future null infinity have recently been resolved and utilized in numerical relativity simulations. With this, gravitational wave models and our understanding of the fundamental nature of general relativity have been vastly improved. In this paper, we review the history and intuition behind memory effects and BMS symmetries, how they manifest in gravitational waves, and how controlling the infinite number of BMS freedoms of numerical relativity simulations can crucially improve the waveform models that are used by gravitational wave detectors. We reiterate the fact that, with memory effects and BMS symmetries, not only can these next-generation numerical waveforms be used to observe never-before-seen physics, but they can also be used to test GR and learn new astrophysical information about our Universe. 

Abstract Type Ia supernovae (SNe Ia), critical for studying cosmic expansion, arise from thermonuclear explosions of white dwarfs, but their precise progenitor pathways remain unclear. Growing evidence supports the “double-degenerate scenario,” where two white dwarfs interact. The absence of nondegenerate companions capable of explaining the observed SN Ia rate, along with observations of hypervelocity white dwarfs, interpreted as surviving companions of such systems, provide compelling evidence for this scenario. Upcoming millihertz gravitational-wave observatories like the Laser Interferometer Space Antenna (LISA) are expected to detect thousands of double-degenerate systems, though the most compact known candidate SN Ia progenitors produce marginally detectable signals. Here, we report observations of ATLAS J1138-5139, a binary white dwarf system with an orbital period of just 28 minutes. Our analysis reveals a 1M_☉carbon–oxygen white dwarf accreting from a high-entropy helium-core white dwarf. Given its mass, the accreting carbon–oxygen white dwarf is poised to trigger a typical-luminosity SN Ia within a few million years, to evolve into a stably transferring AM Canum Venaticorum (or AM CVn) system, or undergo a merger into a massive white dwarf. ATLAS J1138-5139 provides a rare opportunity to calibrate binary evolution models by directly comparing observed orbital parameters and mass-transfer rates closer to merger than any known SN Ia progenitor. Its compact orbit ensures detectability by LISA, demonstrating the potential of millihertz gravitational-wave observatories to reveal a population of SN Ia progenitors on a Galactic scale, paving the way for multimessenger studies offering insights into the origins of these cosmologically significant explosions. 

We present an adaptive-order positivity-preserving conservative finite-difference scheme that allows a high-order solution away from shocks and discontinuities while guaranteeing positivity and robustness at discontinuities. This is achieved by monitoring the relative power in the highest mode of the reconstructed polynomial and reducing the order when the polynomial series no longer converges. Our approach is similar to the multidimensional optimal order detection strategy, but differs in several ways. The approach isa prioriand so does not require retaking a time step. It can also readily be combined with positivity-preserving flux limiters that have gained significant traction in computational astrophysics and numerical relativity. This combination ultimately guarantees a physical solution both during reconstruction and time stepping. We demonstrate the capabilities of the method using a standard suite of very challenging 1d, 2d, and 3d general relativistic magnetohydrodynamics test problems. 

SpECTRE is an open-source code for multi-scale, multi-physics problems in astrophysics and gravitational physics. In the future, we hope that it can be applied to problems across discipline boundaries in fluid dynamics, geoscience, plasma physics, nuclear physics, and engineering. It runs at petascale and is designed for future exascale computers. SpECTRE is being developed in support of our collaborative Simulating eXtreme Spacetimes (SXS) research program into the multi-messenger astrophysics of neutron star mergers, core-collapse supernovae, and gamma-ray bursts. 

Abstract We report the discovery of three ultracompact binary white dwarf systems hosting accretion disks, with orbital periods of 7.95, 8.68, and 13.15 minutes. This significantly augments the population of mass-transferring binaries at the shortest periods, and provides the first evidence that accretors in ultracompacts can be dense enough to host accretion disks even below 10 minutes (where previously only direct-impact accretors were known). In the two shortest-period systems, we measured changes in the orbital periods driven by the combined effect of gravitational-wave emission and mass transfer. We find $\dot{P}$is negative in one case, and positive in the other. This is only the second system measured with a positive $\dot{P}$, and it is the most compact binary known that has survived a period minimum. Using these systems as examples, we show how the measurement of $\dot{P}$is a powerful tool in constraining the physical properties of binaries, e.g., the mass and mass–radius relation of the donor stars. We find that the chirp masses of ultracompact binaries at these periods seem to cluster around ${\mathcal{M}}_c\sim 0.3M_{\odot }$, perhaps suggesting a common origin for these systems or a selection bias in electromagnetic discoveries. Our new systems are among the highest-amplitude known gravitational-wave sources in the millihertz regime, providing an exquisite opportunity for multimessenger study with future space-based observatories such as LISA and TianQin. We discuss how such systems provide fascinating laboratories to study the unique regime where the accretion process is mediated by gravitational waves. 

Abstract We present a discontinuous Galerkin (DG)–finite difference (FD) hybrid scheme that allows high-order shock capturing with the DG method for general relativistic magnetohydrodynamics. The hybrid method is conceptually quite simple. An unlimited DG candidate solution is computed for the next time step. If the candidate solution is inadmissible, the time step is retaken using robust FD methods. Because of its a posteriori nature, the hybrid scheme inherits the best properties of both methods. It is high-order with exponential convergence in smooth regions, while robustly handling discontinuities. We give a detailed description of how we transfer the solution between the DG and FD solvers, and the troubled-cell indicators necessary to robustly handle slow-moving discontinuities and simulate magnetized neutron stars. We demonstrate the efficacy of the proposed method using a suite of standard and very challenging 1D, 2D, and 3D relativistic magnetohydrodynamics test problems. The hybrid scheme is designed from the ground up to efficiently simulate astrophysical problems such as the inspiral, coalescence, and merger of two neutron stars.